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The crystal structure of human chymase complexed with a novel benzimidazole 
inhibitor, TJK002, was determined at 2.8 A resolution. The X-ray crystal- 
lographic study shows that the benzimidazole inhibitor forms a non-covalent 
interaction with the catalytic domain of human chymase. The hydrophobic 
fragment of the inhibitor occupies the SI pocket. The carboxylic acid group of 
the inhibitor forms hydrogen bonds with the imidazole N(e) atom of His57 and/ 
or the 0(y) atom of Serl95 which are members of the catalytic triad. This 
imidazole ring of His57 induces tt-tt stacking to the benzene ring of the 
benzimidazole scaffold as P2 moiety. Fragment molecular orbital calculation of 
the atomic coordinates by X-ray crystallography shows that this imidazole ring 
of His57 could be protonated with the carboxyl group of Aspl02 or hydroxyl 
group of Serl95 and the stacking interaction is stabilized. A new drug design 
strategy is proposed where the stacking to the protonated imidazole of the drug 
target protein with the benzimidazole scaffold inhibitor causes unpredicted 
potent inhibitory activity for some enzymes. 
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1. Introduction 

Chymase is a mast-cell-specific serine protease that is stored 
within secretory granules and released together with heparin 
and histamine in response to allergen challenge or other 
stimuli. Recent studies have shown that chymases possess 
processing activity for biological peptides and cytokines 
implicated in a variety of diseases. For example, the chymases 
of primates and dogs have highly specific angiotensin II 
(Ang II) generating activity, and the results of animal studies 
suggest that chymase contributes to the pathogenesis of 
cardiovascular diseases via Ang II generation (Urata et al, 
1990; Takai et ah, 1997). Based on the results of these studies, 
inhibition of chymase is expected to provide therapeutic 
means for the treatment of these diseases. 

An early approach toward the design of potent inhibitors 
for chymase has been to develop molecules containing acti- 
vated ketones that achieve potency through the formation (or 
likely formation) of covalent adducts with the Serl95 or His57 
residues of the catalytic triad (Aoyama et ai, 2001; Akahoshi 
et ah, 2001). Selectivity in this type of inhibitor becomes 
critical and such inhibitors with large molecular weight tend to 
show a lack of the chances of obtaining oral availability and 
minimal toxicity. This is prompting the search for non-covalent 
reversible inhibitors. There are relatively few reports 
describing inhibitors that specifically and rationally exploit 



non-covalent interactions with the common catalytic residues 
of chymase. 

Recently, we have developed some novel benzimidazole 
derived human chymase inhibitors and determined the crystal 
structures of the human chymase. The benzimidazole inhibitor 
TJK002 (Fig. 1; Yajima et al, 2012) showed potent inhibitory 
activity (K { value 2.24 nM) with respect to human chymase. 
The crystal structure of human chymase with TJK002 was 
determined at 2.8 A resolution. X-ray crystallographic study 
showed that TJK002 forms a non-covalent interaction with the 
catalytic domain of human chymase. The 4-methylbenzothio- 
phen-3-yl fragment of TJK002 occupies the SI pocket. The 
carboxylic acid fragment of TJK002 forms hydrogen bonds 
with the imidazole N(e) atom of His57 and/or the O(y) atom 



OH 




Figure 1 

Benzimidazole derived human chymase inhibitor, TJK002. 
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of Serl95 which is a member of a catalytic triad. This imida- 
zole ring of His57 forms jr-jr-stacking with the benzene ring of 
the benzimidazole scaffold as P2 moiety (Takenaka et al, 
1984). Fragment molecular orbital (FMO) calculation of the 
atomic coordinates by X-ray crystallography showed that this 
imidazole ring of His57 could be protonated with the carboxyl 
group of Aspl21 or the oxyanion group of Serl95 and the 
stacking interaction between the benzimidazole group and 
His57 is stabilized. We propose a new drug design strategy 
using the stacking interaction of the protonated imidazole 
with the inhibitor causing unpredicted potent inhibitory 
activity even for other drug targets. 

2. Materials and methods 

2.1. Crystallization 

The crystallization experiment of human chymase with 
TJK002 was performed using the hanging-drop vapour-diffu- 
sion method. The ligand was added to aliquots of the purified 
protein in a five-fold molar excess. Crystallization conditions 
were similar to those for the PMSF-treated human chymase 
crystal by mixing 1 ul of 6 mg ml -1 protein solution with an 
equal volume of the reservoir solution, which contained 
100 mM sodium citrate (pH 5.5), 15% PEG1500 and 20% 
2-propanol, and equilibrating against 1 ml reservoir solution 
(McGrath et al., 1997). Single crystals grew to suitable 
dimensions in 2-4 d. Prior to flash-freezing in liquid nitrogen, 
the human chymase crystal was transferred to the reservoir 
solution with 30% glycerol and cooled at 79 K. 

2.2. Data collection and structure determination 

Diffraction data were collected on beamline NW12 at the 
Photon Factory (PF) and processed using the HKL2000 soft- 
ware package (Otwinowski & Minor, 1997). Molecular 
replacement was performed using MOLREP (Vagin & 
Teplyakov, 1997) from CCP4 (Collaborative Computational 
Project, Number 4, 1994) with the coordinates of the PMSF- 
treated human chymase (PDB code lklt; the solvent mole- 
cules and PMSF were removed) as the initial model. Refine- 
ment was carried out using the program REFMAC 
(Murshudov et al., 1997). A sample containing a random 5% of 
the total reflections in the data set was excluded for R free 
calculations. After rigid-body refinement, the electron density 
for the TJK002 ligand was clearly constructed using COOT 
(Emsley & Cowtan, 2004). In the final refinement at 2.8 A, 
the crystallographic Rf act0 r an d ^f ree were 26.8 and 32.1%, 
respectively. Statistics of the data collection and final structure 
are summarized in Table 1. Figures were produced using DS 
Visualizer (Accelrys, http://accelrys.co.jp/). 

2.3. FMO and molecular orbital energy calculation 

The atomic coordinates of the complex used for FMO 
calculation used the complemented structure at the defected 
outer loop of the present chymase TJK002 crystal structure 
(PDB code: 4kp0). Construction of the defected loop and H 
atoms were generated using the Discovery Studio software, 



Table 1 

Data collection and refinement statistics for chymase/TJK002. 
Values in parentheses are for the highest-resolution shell. 



Space group P\{l{L 

Unit-cell parameters (A) 

a 56.9 

b 56.9 

c 177.5 

Data collection 

Beamline NW12 (PF) 

Wavelength (A) 1.0000 

Resolution (A) 50.0-2.80 

Total number of reflections 211779 

Unique reflections 27558 

# merge 0.068 (0.422) 

Completeness (%) 100.0 (100.0) 

Multiplicity 15.2 (13.9) 

1/(7(1) 40.7 (5.05) 

Refinement statistics 

Resolution (A) 41.0-2.80 

^factor (%) 24.2 

R fiee (%) 33.0 
RMS deviation from ideal values 

Bond length (A) 0.012 

Bond angle (°) 1.70 



Emerge = E£z'|/z' — where (I) is the mean intensity of N reflections with 

intensities Ii and common indices h, k and /. ^ fact or = S^/||F obs | — k\F C2L \^\lll hk i\F oh ^, where 
F obs and F calc are the observed and calculated structure factors, R fTee is calculated for a 
randomly chosen 5% of reflections and ^ factor is calculated for the remaining 95% of 
reflections. 

version 3.5.0 (Accelrys, San Diego, CA, USA). The 
CHARMM force field implemented in Discovery Studio 
version 3.5.0 was used for the minimization steps. Protein 
structures were optimized with the heavy atoms constrained of 
the crystal structure besides the inserted loop region. Using 
the FMO method, a molecule or molecular cluster is divided 
into M fragments (monomers). Ab initio calculations are then 
performed repeatedly on these monomer fragments in the 
presence of the electrostatic potential created by surrounding 
(M — 1) monomers (Vj), until all the monomer densities 
become consistent. Next, the dimer equations are solved in the 
presence of the electrostatic potential from neighbouring 
(M — 2) monomers (V IJ ). Finally, the total energy of the 
system, E, is written as equation (1) using the total energies of 
the monomer Ej and the dimer E n , 

M M 

E=Y.E U -(M-2)Y.E I . (1) 

The inter-fragment interaction energy (IFIE) in the FMO 
calculations are defined by 

AE IS = {E' I} - % - E'j) + Ti(P u V IJ ), (2) 

where P IJ is the differential density matrix, V IJ is the envir- 
onmental electrostatic potential for the dimer, and E[ and E' u 
are the energies of the monomer and the dimer, respectively, 
in the absence of an environmental electrostatic potential. In 
this study, human chymase is divided into single amino acid 
residues to investigate intermolecular interactions between 
the enzyme and TJK002. The fragmented residues did not 
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exactly correspond to each amino acid residue because 
partitioning of the chemical structure in the FMO scheme was 
performed between the Ca atom and the main-chain carbonyl 
group. Thus, the main-chain carbonyl group of the zth residue 
was assigned to the (/+l)th residue fragment. Using such 
fragmentation, each IFIE is regarded as a residue-residue 
interaction energy. Single-point energy calculations were 
performed by Hartree-Fock (HF) and second-order M0ller- 
Plesset perturbation (MP2) methods (Fedorov & Kitaura, 
2009), using the 6-31G basis set (FMO-HF/6-31G and FMO- 
MP2/6-31G). All FMO calculations were performed with the 
ABINIT-MP (version 5.0) program and the results visualized 
with Biostation Viewer (version 8.0) (CISS Free Software: 
http://www.ciss.iis.u-tokyo.ac.jp/english/dl/). Calculations were 
carried out on a DELL PowerEdge 1950 III (Quad-Core Xeon 
X5470 3.33 GHz x 2) cluster. 

We extracted only TJK002 and His57 coordinates from the 
present crystal structure and calculated molecular orbital 
energy using the quantum chemistry DFT/B3LYP method 
(Jaguar, version7.9: http://www.schrodinger.com/). The basis- 
set function used in the calculation was 6-3 1G**. 

2.4. Determination of K\ 

K { values for human chymase inhibition were determined 
using Suc-Ala-His-Pro-Phe-pNA (Bachem). Enzyme activities 
were evaluated by measuring /7-nitroaniline release from the 
synthetic substrates. K { values were determined by secondary 
plot methods using enzyme activities in some inhibitor and 
substrate concentrations. 

3. Results and discussion 

3.1. Binding conformation and hydrophobic interaction of 
TJK002 

The benzimidazole scaffold inhibitor TJK002 showed 
potent inhibitory activity (K x value of 2.24 nM) with respect to 
human chymase. The crystal structure of the human chymase 
with TJK002 was determined at 2.8 A resolution (Fig. 2). The 
X-ray co-crystal structure of TJK002 bound to the active site 
of human chymase is detailed in Fig. 3. X-ray crystallographic 
study showed that TJK002 forms a non-covalent interaction 
with the catalytic domain of human chymase. In particular, the 
4-methylbenzothiophen-3-yl fragment of TJK002 occupies the 

51 pocket as shown in Fig. 3. It is clear that the potent in- 
hibitory activity that TJK002 showed is derived from this 
effective hydrophobic interaction. The benzene ring of the 
benzimidazole scaffold shallowly occupies the S2 domain and 
could not form an effective hydrophobic interaction to the 

52 site. 

3.2. Oxyanion hole interaction of TJK002 

The various close contacts that lead to important ligand- 
protein interactions are shown in Fig. 4. The alkyl chain of S- 
butyric acid is bent and the carboxylic acid motif of TJK002 
forms hydrogen bonds to the residues which are members 
of a catalytic triad (His57, Serl95). There are known serine 




Figure 2 

Overall structure of the human chymase TJK002 complex. 




Figure 3 

Crystal structure of TJK002 bound to the chymase active site. 

protease inhibitors that possess phosphonate or phosphinate 
groups, which can occupy the active site in the vicinity of the 
catalytic residues Serl95 and His57 (Greco et al., 2007). 
However, phosphonate or phosphinate groups are known as 
undruglike motifs in the viewpoint of the pharmacokinetics 
profile. In our case the druglike carboxylic acid motif of 
TJK002 supplies two O atoms in the binding to the protein, 
which represent a new serine protease inhibitor motif. The 
position of this carboxylic acid motif is occupied by a water 
molecule in the case of another benzimidazole scaffold in- 
hibitor (Lo et al, 2011). 

There are two possibilities from this crystal structure, 
protonated His57 state (Fig. 4a) and neutral His57 state 
(Fig. 4b). In the former state, His57 is protonated by Serl95 
(Kawamura et al., 2011). It is difficult to determine the proton 
position by X-ray analysis. The imidazole ring of His57 forms 
7T-7T-stacking to the benzene ring of the benzimidazole scaf- 
fold as shown in Fig. 5. So far the importance of the S2 site has 
not been noted, but this tt-tt interaction between the benzene 
ring of the benzimidazole scaffold inhibitor at the S2 site 
makes sense because the stacking interaction is also found in 
the complexes crystallized under different pH conditions (Lo 
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Figure 4 

Depiction of the principal interactions: (a) protonated His57; (b) neutral 
His57. 



et al, 2011). The FMO calculation was performed in order to 
understand the features of this 7i-7r-stacking interaction. 



3.3. His57 stabilization with TJK002 by FMO and molecular 
orbital energy calculation 

Fig. 6 shows the IFIE values between human chymase and 
TJK002 in the case of protonated His57 and neutral His57. 
The protonated His57 has a more stable IFIE value 
(-13.48 kcal mol -1 ) than the neutral His57 (+2.94 kcal mol -1 ) 
by binding of TJK002. FMO calculation of the atomic coor- 
dinates by X-ray crystallography showed that this imidazole 
ring of His57 could be protonated through hydrogen bonding 
with the carboxyl group of Aspl02 or hydroxyl group of 
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Figure 5 

7T-7T-stacking interaction between ligand and protein, (a) Perpendicular view; (b) side view. 



Figure 6 

FMO calculation results. Energy in kcal mol -1 . Blue: protonated His57. 
Red: neutral His57. 



Serl95, and the stacking is stabilized in the crystallization 
performed at pH 5.5. The distance shows the strong hydrogen- 
bonding interaction between the Ne atom of the protonated 
His57 and the carbonyl oxygen of the TJK002 (2.6 A, shown in 
Fig. 4). This direct interaction could explain the potent inhi- 
bitory activity of TJK002. In order to understand the inter- 
action energy, a FMO calculation was performed with His57 in 
its protonated and neutral states. The result clearly shows that 
even the neutral His57 shows a slightly repulsive IFIE value 
with TJK002, but TJK002 shows a dominant interaction with 
the protonated His57, similar to Lysl92 of the SI pocket and 
Serl95 as an oxy anion state. This result suggests that the 
protonated His57 is the probable state and the interaction 
with His57 plays an important role in the binding of TJK002 at 
the S2 site at physiological pH in the human body, the same as 
for SI site. A similar stacking interaction between His57 and 
the benzimidazole ring is shown in human chymase complex 
which was crystallized at pH 8.0 (Lo et al, 2011) in spite of the 
carboxyl group interactings with Lysl92 not His57. His57 
should be protonated at physiological pH. 

Fig. 7 shows NHOMO and LUMO orbitals and Table 2 
shows the energy values of NHOMO, HOMO, LUMO and 
NLUMO. LUMO corresponds to the protonated imidazole 
ring and NHOMO corresponds to the benzimidazole ring. The 
energy gap between LUMO and NHOMO is 3.764 eV. This 
might be a possible value for the charge transfer through tt-tt- 
stacking interaction (Wu et al, 2008). The 7r-orbitals of both 
rings have suitable symmetry and 
overlap efficiently. 



4. Conclusion 

We have found a novel, potent and 
non-covalent small-molecule human 
chymase inhibitor. The crystal structure 
of TJK002 bound to chymase shows 

(i) that potent inhibitory activities are 
derived from effective hydrophobic 
interaction to the SI pocket of chymase, 

(ii) drug-like carboxylic acid acts as a 
new serine protease inhibitor motif and 
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Table 2 

Result of molecular orbital energy calculation (DFT/B3LYP method, 
6-31G**). 





Energy (Hartree) 


Energy gap 


NHOMO 


-0.293332 




HOMO 


-0.285330 




LUMO 


-0.155133 


0.138 x 26.2116 = 3.764 eV 


NLUMO 


-0.126900 









Figure 7 

Result of the molecular orbital energy calculation, (a) NHOMO; (b) 
LUMO. 



effectively interacts with the catalytic triad residues, (iii) the 
benzene ring of benzimidazole scaffold forms a 7r-7r-stacking 
interaction with the protonated His57 that could compensate 
for the lack of hydrophobic interaction efficacy between the 
benzene ring of the benzimidazole scaffold and S2 site and 
also be used as the dominant guiding force for the potent 
inhibitory activity of this molecule. We propose a new drug 
design strategy using the stacking interaction of the proto- 
nated imidazole with the inhibitor causing the unpredicted 
potent inhibitory activity even against other drug targets. 
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